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Abstract. We analyze a special class of 1-D quantum walks (QWs) realized using optical multi-ports. We 
assume non-perfect multi-ports showir\g errors in the connectivity, i.e. with a small probability the multi- 
ports can connect not to their nearest neighbor but to another multi-port at a fixed distance - we call this 
a jump. We study two cases of QW with jumps where multiple displacements can emerge at one timestep. 
The first case assumes time-correlated jumps (static disorder). In the second case, we choose the positions 
of jumps randomly in time (dynamic disorder). The probability distributions of position of the QW walker 
in both instances differ significantly: dynamic disorder leads to a Gaussian-like distribution, while for static 
disorder we find two distinct behaviors depending on the parity of jump size. In the case of even-sized 
jumps, the distribution exhibits a three-peak profile aroimd the position of the initial excitation, whereas the 
probability distribution in the odd case follows a Laplace-like discrete distribution modulated by additional 
(exponential) peaks for long times. Finally, our numerical results indicate that by an appropriate mapping a 
universal functional behavior of the variance of the long-time probability distribution can be revealed with 
respect to the scaled average of jump size. 

PACS. 3.67.-a, 05.40.Fb, 02.30.Mv 



1 Introduction 

The quantum walk (QW) is a quantum mechanical model, 
a generalization of a classical random walk. It was in- 
troduced in 1993 [1 2, 3 1 and later found fruitful appli- 
cations as a tool to design efficient quantum algorithms. 
The model of a quantum walk was defined in two dis- 
tinct ways: continuous- and discrete-time. In the former, 
the particles (walkers) are achiral, the Hilbert space is 
spanned by the discrete position states corresponding to 
vertices of a graph. In the discrete-time case, the introduc- 
tion of chirality is unavoidable. The Hilbert space corre- 
sponding to the chirality has its dimension equal to the 
number of possible directions of a step. 

A possible experimental implementation of a classical 
random walk is the Gallon board (also known as Quin- 
cunx). Here a large number of balls (walkers) fall through 
the board, changing their direction randomly on period- 
ically arranged pins and forming so a binomial distribu- 
tion of their final position. A quantum analogy of the Gal- 
ton board (and one possible implementation of the QW) 
is shown in Fig. [T] where the walker is a coherent light 
pulse moving through a medium with periodical bound- 
aries that split the signal; finally, there are detectors at the 



end which represent the quantum equivalent of the bins 
in the classical model. 

The spectrum of investigation of the QW (and its 
modifications) is broad. The original idea was presented 
in fTT'S] and since then a few review papers have been 
published [4,5 6 1. Recently, there have been studies of 
QW allowing Levy noise in the model. The latter is in- 
troduced via randomly performed measurements with 
waiting times following a Levy distribution JZHSl- The 
most relevant paper for the present work is 1 9] in which 
the properties of the one-dimensional continuous-time 
QW in a medium with static and dynamic disorder are 
examined. Other studies of the QW have focused on 
the meeting problem of two particles 1 10 11 J. Recurrence 
properties of the walker have been investigated in II11II12I 
|13|. Moreover, localization of the walker has been stud- 
ied in [14]. Finally, a theoretical investigation of the QW 
in random environment has been performed in [15 16|, 
where a robust mathematical definition of a random en- 
vironment is provided. 

The simplest analytic task in the study of quantum 
walks is to determine the functional expression of the 
probability to find a particle at a certain location at a 
certain time. General analytic calculations of quantum 
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walks are not known. However, asymptotic solutions of 
several quantum walk models have been found using 
path integral techniques in [14,15 ,16 ,17,18 ,19 ,20 ,21,22] 
and in 1 10 1 1 1 using Fourier transform. 

Since their introduction, many experimental groups 
have tried to implement quantum walks. A number of 
successful realizations of one-dimensional QW have been 
reported in optical lattices |23|, trapped ions |24J, and 
cavity QED |25|. More recently, additional realizations of 
a quantum walk using atoms in an optical lattice [26 1, 
trapped ions [27.28 1 and photons [29| have been an- 
nounced. 

Recent studies of quantum systems II30II31I with ran- 
dom potentials have shown that localization of the parti- 
cle can occur as a result of the randomness. Focusing on 
quantum walks, Ribeiro et al. |32| have used two different 
corn operators switched according to the Fibonacci series 
and they have observed localization in the system. Yin 
et al. [9| have numerically simulated the continuous-time 
QW on a line and they have observed Anderson local- 
ization only in the case of static disorder, while d5mamic 
disorder leads to decoherence and a Gaussian position 
distribution. In addition, effects of spatial errors have 
been studied by Leung et al. in [33 1. They have examined 
QWs in 1 and 2 dimensions on networks with percolation 
where the missing edges or vertices absorb the walker, 
leading to topological randomness of the graph. 

The simplest discrete QW is described by the action of 
the coin and the step operator. Much attention was paid 
to the alternation of the coin operator - position or time 
dependent coin and its implications on the walker dy- 
namics have been extensively discussed. However, little 
focus was given to changes of the step operator. Our anal- 
ysis takes a step in this direction. While assuming a con- 
stant coin we study changes in the step operator. When 
assuming a Galton board realization this amounts to a 
change in the connectivity between the layers of beam 
splitters forming the walks. 

In the present paper, we focus on the QWs where sig- 
nals can jump to a distant location, which is a general- 
ization motivated by Levy flights in classical mechanics. 
The jump may be caused by an rnhomogeneity of the 
material, spatial proximity of non-neighbor channels or 
scrambling in the topology of the network (for instance 
relabeling of input-output label of multiports forming the 
network realising the QW). 

We study two basic modifications of the QW. In the 
first case, which we call dynamic disorder, the jumps are 
prepared as independent and identically distributed in 
time, whereas in the second case, called static disorder, 
the positions of the jumps are perfectly correlated. We in- 
vestigate the problem using computer simulations em- 
ploying the Zarja library [34| Q and its offspring library 
focused on the QW0. 

The structure of the article is as follows. In Section 2, 
we motivate our study from an experimental point of 




Fig. 1. A schematic of an interferometer simulating material and 
implementing a random walk. The red lines are the connections 
of multi-ports along which the signal is transmitted. The light 
blue half-circles represent detectors. The blue arrows are two 
selected input channels. Every horizontal layer of the multi- 
ports is accessed simultaneously. 

view. Then we define a model with next neighbor inter- 
actions only. In Section 3, we show the results of the sim- 
ulations based on Monte Carlo method. In Section 4, we 
discuss our results and draw conclusions. Finally, we de- 
scribe in the appendix the set of operators used in the 
simulations and an algorithm to sample from the appro- 
priate probability distribution. 

2 Definition of a Quantum Wall< 

We model a sequence of optical layers by an array of 
multi-ports, or beam splitters (Fig. |2), forming so a large 
interferometer, see Fig. [1] - this is a quantum analog to 
Galton's board (Quincunx). Due to the regular structure 
of the interferometer, it is natural to treat the temporal 
evolution of the system in discrete steps, discretized by 
the time needed for an excitation to travel the distance 
between two consecutive layers at a constant angle. If 
we let this angle approach the right angle, we obtain so- 
called static disorder (see Fig.|3}i,e), where the same set of 
beam splitters is used repeatedly in every time step. An 
excitation enters the system of multi-ports at one selected 
position and spreads due to the coherent interaction with 
the beam splitters. A number of detectors is placed in the 
interferometer so that the signal hits one detector in every 
possible path after a given number of interactions with 
the medium. The resulting scheme is an implementation 
of the QW. 

In the following, we introduce two ways to describe 
the system. First we introduce the standard notation for a 
discrete-time QW. We define the basis states of the system 
as states localized at the beam splitters just before or after 
the interaction takes place. Therefore, every basis state is 
described by a position within the beam splitter layer (an 
integer number) and a chirality, which can take values L 
or R for left and right, respectively. Formally, we define 
the state space Ti of the system as 

n = Us®nc (1) 



^ |http: //sourcef orge .net/pro jects/zar ja/| where a "position" space V-s and a "coin" space V-c 

^ |http: / /sourcef orge . net/pro jects/ quant umwalk/| (the name stemming from the idea of a random walker 
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Fig. 2. Action of a beam splitter on the input states a'{t) and 
V^{t) which are transformed to the output states c'{t + 1) and 
d'*' + 1) using the matrix M. 



tossing a coin to decide the direction of his next step) are 
defined as 

(2) 



T-Ls = Span{|n) | « e Z} 



and 



?^C=Span{|L),|R)}. 



(3) 



We define all the numbered states | n) , as well as the states 
|L) and |R), to be normalized and mutually orthogonal. 
Therefore, Hs is isomorphic to the space of quadrat- 
ically integrable complex sequences and "He to C^. The 
basis states are then constructed as tensor product basis 
states \ n,R) and \n,L), where « e Z. 

Following Konno |16|, we define that the state of the 
system at time Hn a random environment co is described 
by the positive semi-definite density matrix Pu){t) on 
% having Tr pco{t) = 1. Thus the density matrix of the 
system p{t) at time f is 



p{t)= / pco{t)-F{du;), 



(4) 



where symbols O and P will be specified later in section 
2.2. 

The evolution of the system by one timestep from 
t to t + 1 for one particular setup of environment co is 
described by a unitary operator Ua>{t t + 1) which acts 
as follows: 

Pco{t + l) = Ua;{t^t + l)-pa,it)-Ulit^t + l) (5) 

just the same way as typical propagator on density ma- 
trix in quantum mechanics in, e.g., 1351 . 



2.1 Quantum Walk 

The QW described, e.g., in ||4| assumes that the condi- 
tions of the environment are stable and thus it holds 
p{t) =pco{i) for all CO e O. If we set up the initial den- 
sity matrix p[Q) as a pure state ji/'o) G Ti using p(0) 

('/'ol/ then the system is in a pure state \ipj 
\jQ^ {i t + 1) |i/'o) (the product expanded in the proper 
time ordering) at every timestep T and the appropriate 
density matrix is p{T) = \ipj){^j\. Moreover, the original 
QW was defined with fixed unitary operator U'^^ {t — > 



ttt=0 



f + 1 ) in the form of a composition of two unitary opera- 
tions as follows: 



(6) 



Thus, the evolution of the walker at time T is Ii/'t) = 

(UQ^) |!/?o>- The operator C is called a coin operator 
and describes the transformation induced by a beam 
splitter. Here the position of a localized state stays un- 
changed and the operation acts only on the coin state as 
C = I ® M, where M is a unitary operation transform- 
ing the probability amplitudes due to a partial reflection 
on the beam splitter. We assume for simplicity that all the 
beam splitters have the same physical properties and per- 
form a Hadamard transform on the input states. 



M 



(7) 



The operator S represents the propagation of an excita- 
tion in the free space between the beam splitters. Hence, 
the chirality of a basis state does not change but the posi- 
tion is shifted by ±1, depending on the coin state. We can 
express S as 



S= ^(|n + l,R)(n,L| 



l,V){n,R\ 



(8) 



It is easy to show that the sum converges and defines a 
unitary operator defined on all the state space T-L . 

If the initial state of the walker is one of the basis 
states, for example |0,R), it evolves under the operation 
U so that in terms of a complete measurement in the 
position space, the probability spreads to both sides from 
the starting position. However, it is bounded between the 
positions —t and f as it can't change by more than 1 in 
either direction in any step. Moreover, due to the fact that 
a transition by 1 has to be done in every time step, the 
walk is restricted at each time f to a subspace spanned 
by the basis states for which the position shares the same 
parity with f . 

Alternatively, we can describe the path of the walk 
in terms of the edges of the underlying graph instead of 
its vertices, following the physical trajectories of the exci- 
tation and eliminating the need to describe the propaga- 
tion between beam-splitters. In the following, we will call 
these edges "channels". Due to the fact that every chan- 
nel lies in between two consecutive positions of the beam 
splitters, as well as to avoid confusion in the notation, 
we will denote the channels by half -integer numbers. A 
notation based on the channel formalism can be intro- 
duced and mapped to the previously defined state space 
in more ways. One such possibility is to identify the chan- 
nel as a given time t with the state in which the excitation 
is at the end of the propagation, just before hitting an- 
other beam splitter (or a detector). This time-dependent 
mapping is given by the formulas 



\n+\;t) 




if t is even, 
if t is odd 



(9) 
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for even n and 



\n + i;t) 




if t is even, 
if t is odd 



(10) 



for odd n. One can verify that such a mapping respects 
the parity rule and covers all the subspace of H that is 
actually used by the quantum walk if we assume that the 
initial state was \0,R) or |0, L). In particular, these two 
initial states are denoted | — j;0) and \j;0) in the channel 
notation, respectively; cf. Fig. [1] for a visualization. Since 
the instantaneous chirality of an excitation is uniquely 
given by the position of the channel and the time, we 
do not need to specify the coin degree of freedom ex- 
plicitly in this approach. The "corn toss" and "step" are 
merged into one unitary operation which mixes neigh- 
boring channels in pairs. 

The above discussion gives two equivalent ways to 
describe QW on a line. Throughout this work, we will 
use the latter approach as it makes it much simpler to 
describe the jumps in the network. 



combinations of jump operators P^, exchanging signals 
in channels /, and + /, having the matrix form 



/lO . 
' 1 . 


oo 


oo 


. 0\ 
.00 


6 . 


.' . 


. 1 . 


. 6 


6 . 


. 1 . 


.' . 


. 6 


.00. 

\o . 


oo 


oo 


.' 1 6 , 

. 1/ 



(12) 



i.e., that of a transposition operator. = 2^ is the cr-field 
defined on 



n = {E,p,,p,p,,...} 



(13) 



and ¥ : J- ^ [0, 1] is the probability measure, specified by 
the probability of elementary events, defined as 



(14) 



2.2 Quantum Walk with jumps 

We will assume that the quantum walk is disturbed by 
random topological errors, see Fig.|3j which are modeled 
by random changes of connectivities between the multi- 
ports forming the network in Fig.[TJ We focus our analysis 
on two distinct situations. First we will assume that we 
deal with the repetition of random but stationary errors. 
Stationary errors mean that in each layer the same jumps 
appear. One could also realize static disorder with a sin- 
gle set of beam splitters, if the incident beam is parallel 
to the layer of beam splitters and repeatedly sent through 
them (see Fig.|3]i,e). The other situation refers to the case 
when in each layer jumps are generated independently 
We refer to the latter situation as dynamic disorder (see 
Fig.|3i,b,c). 

In the QW with jumps we assume that the unitary 

operators Ua,{t ->t + l) = ut"'^{t ^ f + 1) in Eq.|5]are 
random (the set of such operators forming a probability 

space of random unitary operators) and that Ul^^"'^{t 
t + 1) can be written in form of a joint action of two 
unitary operations. 



ur'{t -^t+i) = sr'it) ■ uS'^it ^ f + 1), (11) 

where {t -^t + 1) is the evolution operator of the ID 
QW in a clean media without errors, defined by Eq. |6] 

We define that the set of all unitary operators S^",'^^ (t) 
forms the probability space S(n,J-',P). For the sake of 
simplicity of the following discussion and the numerical 
simulations, we make the set Cl finite by replacing the 
infinite walking space by a cycle of size N with a periodic 
boundary condition. If N is chosen sufficiently large, 
this imposes no restriction on the validity of the results. 
The sample space O is then the set of all the possible 



We define tr{n) as the number of transpositions of in- 
dexes forming permutation tt, N is equal to the size of 
the system and consequently to the dimension of jump 
operators Pj., p is probability that one pair of errors with 

distance / occurs and finally Z=(l + ( — p)sj is nor- 



malization where g = gcd(N,7') (gcd stands for greatest 
common denominator). 



3 Simulation of tlie model 



The disorder of the system is expressed by the set of uni- 
tary operators displacing the signal to a distant location. 
In general, the system can undergo either static or dy- 
namic disorder. The system with static disorder is prop- 
agated by random unitary operators fixed during the re- 
alization, i.e., the operators are correlated in time. On the 
other hand, the evolution of the system undergoing dy- 
namic disorder is governed by the operators that are in- 
dependently and identically distributed in time. 



3.1 Static disorder 



The random unitary operators propagating the system 
with typical structure of jumps according to Eq. [12] can 
have due to mapping |9] and [TO] one of two fundamental 
distinguished forms, according to the parity of the length 
of the jumps induces by the errors in the media. Jumps 
of even lengths j do not change chiral state of the walker 
in Hilbert state H in contrast to jumps of odd lengths / 
which swap chirality of the walker from state L to R and 
vice versa. 
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Fig. 3. Schematic picture of different jumps of the signal. 
Scheme (a) shows a part of interferometer without jumps. Two 
schemes in the middle, namely (b) and (c), show dynamic dis- 
order pattern where (b) and (c) shows even (; = 2) and odd 
= 3) size of jump, respectively, which can emerge indepen- 
dently at every level. Scheme (d) shows an interferometer with- 
out jumps for static disorder and a jump of size ; = 2 is illus- 
trated in scheme (e). 

3.1.1 Odd jumps 

The only tunable parameter that affects the evolution of 
the walker is the probability p that one error (jump with 
distance /) occurs. The value of the parameter p close to 
should produce typical chiral ID QW pattern, which 
is clearly visible on the top of Fig. S) Here the QW was 
initiated in a localized state |0 + j;0) of the walker, caus- 
ing an asymmetrically distributed walk. Higher values of 
p produce totally different pattern where high-frequency 
oscillations of probability distribution of positions of the 
walker are suppressed. In the middle and at the bot- 
tom of Fig. 3] the probability distribution of the positions 
of the walker shows Laplace distribution modulated by 
Laplace distributed peaks with distance / between neigh- 
boring peaks (clearly observable as the small triangular 
peaks modulated on a bigger structure in Fig.|4). 

Previous observations of the fundamental change in- 
duced by the variation of p is supported by Fig. |5] show- 
ing evolution of the probability distribution. The top part 
shows the typical structure of evolution of the probabil- 
ity distribution of the walker — the quantum carpet |36 j — 
with additional quantum carpets on the border of the 



Probability distribution among channels 




-200 -100 100 200 



Channel 

Fig. 4. Typical patterns of probability distribution of position 
of the walker among channels after T = 200 steps for ; = 21, 
from the top to the bottom p = 0.01, p = 0.05 and p = 0.2. On 
the top, the typical structure of the QW without jumps that is 
perturbed by formation of an extra central peak. In the middle, 
we see formation of Laplace-like tail and 3 peaks with distance 
j. Central peak is located at the position of the initial excitation 
of the system at T = 0. On the bottom, the whole distribution 
follows an overall exponential decay which is modulated by 
exponential peaks with a distance j = 21 between maxima. 
The graphs were obtained using Monte Carlo method after 
R = 20000 runs of the randomized evolution. 



main one which were formed at initial stages of evolu- 
tion. For small values of p, the interference is not strong 
enough to change the pattern of distribution. In contrast, 
in the second case displayed at the bottom of Fig. |5l 
where p = 0.5, we observe a typical structure of equidis- 
tant peaks separated by valleys of width j. The peaks are 
formed early during evolution and they do not change 
their positions later. 

To support the idea that the walker froze in the close 
neighborhood of a few preferred locations we analyzed 
the evolution of entropy and variance of the probabil- 
ity distribution of the position of the walker. The evolu- 
tion of the variance, shown in Fig. |6j visualizes the ob- 
servation from the previous paragraph. For a probability 
p close to (unperturbed QW), we see a clear ballistic 
diffusion of the walker (tJ^ ^ t-^) — in contrast to this we 
stress that for a classical random walker cr-^ ^ t. Increas- 
ing p, we observe sub-ballistic diffusion and finally, in the 
range close to 1/2, the variance tends to finite constant. 

Next, we turn our attention to the evolution of the 
classical (Shannon) entropy of the probability distribu- 
tion of the positions of the walker, which is shown in 
Fig. [71 The evolution of the classical entropy follows an 
analogous behaviour as described in the case of variance 
in the previous paragraph. Thus, we can observe for the 
probability p close to unperturbed QW an increase of 
classical entropy proportional to In(f). However, rising 
probability p causes slower and slower increase in en- 
tropy in time and for the case p = 0.4 we observe a satu- 
ration of entropy. 
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Fig. 5. Probability distribution among channels for T = 200, 
R = 20000 and ; = 11; from the top to the bottom p = 0.01 and 
p = 0.5. On the top, the maxima of probability move outward 
similarly to the standard case of the QW. On the bottom, the 
maxima are formed and they stay at the point of original for- 
mation. 

The conclusion of the above paragraphs on odd-sized 
jumps in the QW is that we clearly observe localization 
of quantum walker in static disorder media described 
by the probability space of random unitary operators 
S{ri,J^,P) parametrized by probability p. Moreover, the 
numerical results suggest that taking separately odd and 
even timesteps the probability distribution of the posi- 
tion of the walker converges to a stationary distribution: 

{n + ^;2t\pi2t)\n+ P;-" [n + 1) (15) 

(n + l;2f + l|p(2f + l)|n + l;2f + 1) 

Pp^™(n + 2) and V°^^(n + j) are universal distributions 
for odd and even timesteps for probabilities p in range 



Evolution of vaiiance 




1 3 10 30 100 200 

Time 

Fig. 6. Evolution of variance of the system f or T = 200, R = 
20000 and / = 11. Transition from the case of ballistic diffusion 
of the signal for probabilities p close to to the case of sub- 
ballistic diffusion, demonstrated for p = 0.05 can be observed, 
and finally we observe for p = 0.4 that variance tends to a 
constant, i.e., diffusion coefficient approaches and the walker 
ceases to "spread" — a characteristic property of localization 
(37). 



Evolution of entropy 




1 3 10 30 100 200 

Time 



Fig. 7. Evolution of classical (Shannon) entropy of probability 
distribution of positions of walker for T = 200, R = 20000 and 
/ = 11. Classical entropy is a rising function of time t for p = 
0.01 where the fimctional form is proportional to Int. On the 
other hand, increasing probability p causes "freezing" of the 
evolution, as can be observed for p = 0.4. 

close to 2 and the size of the region of convergence is 
also dependent on /. 

Let us consider the fundamental properties of the 

asymptotic probability distribution P^J" (n + j^; vari- 
ables of the probability distribution P^'^^'' + 2) O'^'^ 

timesteps are different but general properties are shared. 
We concluded from Fig. |4] that we observe a formation 
of an overall Laplace distribution modulated by Laplace 
distribution of peaks, both in the form 

P(x) = Cexp(^-ti4^, (17) 

where pi is the mean value of the distribution and a is 
related to variance via Var F{x) = 2- a^. Our aim is to 
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Dependence of modulated peak on probability 



Dependence of modulating peak on probability 
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Fig. 8. Fit of the irwerse parameter | of the probability distri- 
bution of position of the walker for T = 200 and R = 15000 
with various odd sizes of jump j taking x-axis as x = p ■ j and 
y = 1/a. We observe the formation of a U-shape function for 
constant sizes of jump / and changing x = p ■ j with a minimum 




at Xfi 



= 2. 



0.2 0.3 
Probability 

Fig. 9. Fit of the inverse parameter ^ of the central peak of the 
probability distribution of the position of the walker for T = 200 
and R = 15000 with various odd sizes of jump j. Increasing 
probability causes linear increase of the inverse parameter ^ 
with p when / is fixed. Moreover we observe a universal de- 
pendence on p, independent of the length of jump j. 



estimate the inverse parameter i of Laplace distribution 
in two cases for 

- the whole distribution (with mean at the point of 
injection), 

- the modulated peaks (with mean at the center of the 
peak). 

We focus on the shape of the whole probability distri- 
bution. The plots in Fig. |8] suggest a U-shape function of 
the fitted inverse parameter 1 / a of the Laplace distribu- 
tion for all parameter values. We vary the probability p 
and we connect points for the same jump sizes /. Defin- 
ing the X-axis as x = p ■ j we put all minima at a fixed po- 
sition. Thus, for the minima p ■ j = x^/n holds where an 
approximation of the constant is Xm,„ = 2. The values of 
minima of the fitted inverse parameter 1/a, reached for 
p = x„ij„/ j, form an increasing function of jump size /. 
The size of the central peak is /, 7/2 on the left and on the 
right from the maximum at 0. Estimation of the inverse 
parameter l/o of the central peak is shown in Fig.|9l The 
shape of the central peak is an increasing linear function 
of probability of jump p independent of the size of jump 
j. Thus the central peak becomes steeper and steeper with 
increasing p — localization of the walker becomes more 
evident. This is true not only for the central peak but it 
holds true for the other peaks as well, see Fig. S] 

Let us look at the dependence of entropy of probabil- 
ity distribution p (j^ ~^ to fi^d ^ walker at position 

n + J. We plot it in Fig. [10] where the entropy was mea- 
sured by both extensive and non-extensive measures. In 
the first case, the extensive measure is classical (Shan- 
non) entropy. The non-extensive one is the (^-entropy in- 
troduced by Tsallis, which for a particular cj = 1 reduces 
to the classical entropy (for more about (j-statistics, see 
II38II39||401|4TB2| ). Both entropies are decreasing functions 
for increasing disorder measured by p. This leads to a 



counter-intuitive statement that increasing classical dis- 
order organizes quantum system even if measured by a 
non-extensive entropy (in classical random walk longer 
jump causes increasing variance and increasing entropy 
as well). Moreover, the (j-entropy of the probability distri- 
bution of position of the walker for the parameter value 
c] = 2 brings all curves corresponding to different sizes of 
jump j on one single curve. This means that if we take 
into account non-additivity of the system quantified by 
the parameter 1 1 — | (taken form the theory of the non- 
additive (j-entropy), we can map the QWs between each 
other for different sizes of error j holding probability of 
occurrence of pair of error positions p constant for q = 2. 

Finally, the last studied variable for the QW with 
static disorder was the variance of position displace- 
ments Var ^p,j, plotted in Fig. [11] on the left and with 
rescaled axes on right-hand side. The variance forms a 
U-shaped function of p with moving minima for differ- 
ent but fixed /. Rescaling the axis x = p ■ f where a = 1.04 
and axis y = ■ Var Pp y where /3 = 1.67, we clearly see 
from Fig. [11] right inset, that there is a universal U-shape 
function fl^^{^) fulfilling 

/o*ddW=r^-^«''lPx,r (18) 

Due to relation a ~ 1 and the U-shape of function f*^^ {x) 
we can conclude that the overall variance of probabil- 
ity distribution of the walker's position is strongly cor- 
related with the fit of the inverse parameter 1/a (of the 
Laplace distribution) of whole probability distribution. 



3.1 .2 Even jumps 

Let us now assume that the ensemble of the random uni- 
tary matrices is still parametrized by the probability p 
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Entropy of the system 



5 10" 



Dependence ot variance on error 

,3 





0.1 0.2 0.3 0.4 0.5 




0.1 0.2 0.3 0.4 0.5 

Probability of error 

Fig. 10. Classical entropy of the probability distribution is on the 
top, Tsallis entropy for q = 2on the bottom for parameters set to 
T = 200 and R = 15000. Classical entropy, on the top, decreases 
for increasing p but the shape of the dependence is different 
for different values of j. The shape of the decrease becomes the 
same for Tsallis entropy with parameter q = 2 as seen on the 
bottom. 

that the pair of erroneous positions occur, but the dis- 
tance of a pair of errors j is an even number and due to 
the mappings |9] and [10] the chirality of walker does not 
change. The typical formation of Laplace-like distribu- 
tion as in the case of odd jumps is not present, but in- 
stead we observe a 3-peaked structure, as seen in Fig. [121 
(left-right asymmetric due to the initial condition) mod- 
ulated by a periodic function that has its period equal to 
the length of jump To emphasize the difference between 
odd and even we conclude that both cases form a lo- 
cated peak at the initial position of the walker. However, 
odd sized jumps ; cause Laplace-like overall distribution 
while the case of even can form a 3-peaked distribution 
where the central (localized) peak is sharp but the two 
others are broad and peaks with distance / are modulated 
on the overall structure. 

3.2 Dynamic disorder 

The static disorder analyzed in the preceding sections as- 
sumes random, but time-correlated emergence of pairs 




I 1 ' 5 I ■ • • ■ 1 

0.1 0.2 0.3 1 2 3 4 5 6 

Probability of error x 



Fig. 11. Dependence of variance of the probability distribution 
for odd j on p for T = 200 and R = 15000 is shown. The typical 
U-shape function can be observed with shifted position of the 
minimum. On the right, there are the same data with reseated 
X-axis where x = p ■ j" where a = 1.04 and y-axis where y = 
i^PVar{i) where /3 = 1.67. The data fit on the same U-shape 
function reaching their minimum at the same place. 



Probability distribution among channels 




-400 -200 200 400 



Fig. 12. Probability distribution among channels for T = 200 
and R = 15000 in the case of even sized jumps (j = 20). The 
probability distribution forms a 3-peak structure that is modu- 
lated by additional periodical small peaks which are separated 
by;. 



of errors that are expressed in the set of random uni- 
tary operators S(n,.F,P). In contrast, dynamic disorder 
assumes independent and identically distributed ran- 
dom unitary operators from the same probability space 

s(n,j',p). 

Let us consider on the probability distribution of the 
position of the walker, plotted in Fig. [13] on the top. We 
observe decoherence of the quantum walker forming a 
distribution reminiscent of a Gaussian distribution mod- 
ulated by ID QW patterns for small probabilities of jump 
p. On the other hand, large p causes a modulation by val- 
leys with distance / between peaks. The functional de- 
pendence of the standard deviation of the probability dis- 
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Distribution of probability among channels 
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Dependence of standard deviation on error 




0.05 0.1 0.15 0.2 



Probability of error 

Fig. 13. On top, probability distribution among the channels 
for dynamic disorder for T = 100, R = 5000 and j = 40. Prob- 
ability distribution can be fitted to a normal distribution where 
residual patterns of the QW are present for p = 0.01; for p = 0.2, 
formation of periodical valleys modulated on overall distribu- 
tion (similar to static disorder with an odd size of jump) can be 
observed. On the bottom, functional dependence of standard 
deviation of the probability distribution of the position of the 
walker on probability of error p and its linear fit. 

tribution of position of the walker on probability of jump 
p in Fig. [13] on the bottom shows a linear behavior. 



4 Discussion and Conclusions 

We have defined a modification of the QW in ID where 
the environment causes long but fixed-size jumps that 
emerge with a constant probability — the model is no 
longer deterministic but stochastic, depending on a ran- 
dom variable. In one step, first a unitary QW coin and 
shift operator act on the state, then a stochastic dis- 
placement operator generates jumps. We have studied 
two classes of QW with disordered connections between 
beam-splitters. First, d5mamic disorder model shows de- 
coherence that leads to Gaussian distribution modulated 
by residual patterns of QW (for small probabilities of 
jump p) or by valleys (for large probabilities of jump p). 
The standard deviation of the position of the walker is a 



linear fimction of p. In the second case, the QW is per- 
turbed by static disorder. The behaviour of the model de- 
pends on the size of jump and we have investigated odd 
and even jump sizes separately. Even jump sizes cause 
localization of the walker at the initial position and two 
other broad peaks modulated by oscillations with a pe- 
riod of the size of jump /. The focus of the paper, however, 
was mainly given to jumps with odd size. In this case, the 
evolution of the variance shows a transition from the bal- 
listic diffusion to a no-diffusion regime with increasing p 
and this observation is supported by numerical calcula- 
tion of classical entropy which changes from logarithmi- 
cally increasing regime to a no-growth regime. 

In addition, the probability distribution of positions 
of the walker changes from a pattern typical to the QW 
to Laplace distribution modulated by Laplace distributed 
peaks separated by the size of jump The formation of 
Laplace distribution depends on the probability p that a 
pair of errors occurs. Moreover, the model exhibits the 
unusual property that classical disorder in quantum sys- 
tem can decrease, i.e., we observe a decrease of both the 
classical (Shannon) and (j-entropy (with q = 2) of mea- 
surements. Finally, our numerical calculation shows that 
using non-additive t^-entropy with q = 2 there is an uni- 
versal functional dependence of (j-entropy on p for ar- 
bitrary j. The next part of out investigation was turned 
to the variance of the probability distribution. Our nu- 
merical results indicate that there is an universal func- 
tional form of variance of the probability distribution of 
the position of the walker — ^variance multiplied by j^^ 
is an universal U-shaped function of one variable p ■ j"^ 
where a = 1.04 and /3 = 1.67. The functional dependence 
of the universal function shows a minimum which sepa- 
rates two regimes — one with decreasing variance and the 
second with increasing variance. To put this result in the 
broader context of Complex Systems and Game Theory, 
we note that similar behavior of variance of attendance 
has been observed in Minority Game, exhibiting dynam- 
ical phase transition, see, e.g., 1341431 . 
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A Properties of the set of operators 

Let the unperturbed walking space be represented by a 
cycle graph with N vertices where every vertex repre- 
sents a channel in the main text. The errors in the network 
are represented by swapping the walker's probability 
amplitudes between vertices labelled i and i + / mod N 
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for a fixed ;', which happens with a relative probability 
p. With a relative probability 1 — p, a vertex is left intact. 
Finally, a vertex already exchanged with another one can- 
not be used for another transposition in the same permu- 
tation. 

These conditions give the probability of a permuta- 
tion TT in the form 



P(7r) 



1 



Mn) 



(1-p) 



N~2-tr(n) 



(19) 



where tr{n) denotes the unique number of independent 
transpositions forming the permutation n. A factor of 2 in 
the exponent of 1 — p is present due to the fact that every 
transposition reduces the number of unused vertices by 
2. Finally the factor Z^^' is the normalization constant. 
The constant Z^^' is computed as 



L-J 

zw=f:N,p'^(i-p) 

k=0 



N-2-k 



(20) 



where N/^ is the count of all possible permutations formed 
by exactly k non-incident transpositions of size /. 

Lemma 1 Let N >2, let j and N are relatively prime. Then 

zW = l + {-pf. (21) 

Proof Due to the relative primality of j and N, we can 
relabel the vertices by indices through N such that ver- 
tices with successive indices have a distance of / in the 
original numbering. This way, we can reduce the prob- 
lem of finding N/f to a combinatorial problem of finding 
the number of fc-element subsets AC {0,1,...,N — 1} sat- 
isfying the following conditions: 



(a) for all < f < N 

(b) {0,N-1} ^ A. 



In order to find Nj^, we discuss two disjoint cases: 

(i) Let us count the subsets which do not contain N — 1 
as an element. For each such set A = {ai,a2,a^, 

%}/ fli < ^2 < ■ ■ ■ < we denote A = {fli,fl2 — 
l,fl3 — 2, .. — A: + 1}. This is a one-to-one map- 
ping, reducing the problem to finding fc-element 
subsets of N — k elements without any additional 

restriction. This gives {^^'^) possible subsets. 

(ii) Let now N — l£ A. Then the other k — 1 elements of 
A must lie in {1,2, . . . , N — 2}, with no two of them 
successive and N — 2 excluded. This is a variant of 
the subproblem (i), giving (^^_^]~^) possibilities for 
A. 



Adding these results, we find that 



N-k 
k 



N-k-1 
k-1 



(22) 



In order to calculate Z^^^^, we find the generating 
function 



+ 00 



F{x) = ^ z(N)x~ 



N=0 
+00 Lt 



- E E^kXpH-^-pf-'V. (23) 

N=2 A:=0 

Here, for simplicity, we generalize |22] also for N < 2 
and we define by convention (~J) = 0. Using standard 
methods, we obtain for the sum 



1 + px^ 



X + xp — px 



2- 



(24) 



This result can be decomposed into partial fractions as 



X 1 + px 



1, 



(25) 



from which the power series can be derived quickly as 



+ CO +CO 

Fix)=Y:x^+Y:i-pxf-i 

N=0 N=0 



-1+E i + i^"- 

N=Ov^^ ^ 



(26) 



QED. 

If N = 2, the formula [21] cannot be used. Indeed, 
computing Z^^' manually gives 



Z(2) = p+(l-p)2 = l_p. 



(27) 



This is because |22] gives an incorrect result for A: = 1 in 
this case. In practical situations, however, N ^2. 

Lemma 2 Let N and j be positive integers such that j < N, 
g = gcd(N,7) and f > 2. Then 



zw=(^i+(-p)fy 



(28) 



Proof If ^ = 1, then j and N are relatively prime and we 
can use Lemma[T]to obtain the result. In the case g > 1, we 
first split the set {0, 1, 2, . . . , N — 1 } into g modular classes 
(mod g). According to the definition, the selection of 
the errors can be done independently on each of these 



subsets. Therefore, as y is relatively prime to we can 

use Lemma [T] for each of these classes and multiply the 
partial results to obtain the total partition fimction in the 
form stated by the Lemma. 



H. Lavicka et al.: Quantum Walk with Jumps 



11 



References 

1. L. Davidovich, N. Zagury, and Y. Aharonov. Phys. Rev. A, 
48:1687-1690, 1993. 

2. S. Gudder. Quantum probability. Academic Press Inc., CA, 
USA, 1988. 

3. G. Grossing and A. Zeilinger. Complex systems, (2):197-208, 
1988. 

4. J. Kempe. Contemporary Physics, 44:307, 2003. 

5. D. Meyer. /. Stat. Phys., 85:551-574, 1996. 

6. S.E. Vanegas-Andraca. Quantum Walks for Computer Scien- 
tists. Morgan and Claypool, 2008. 

7. A. Romanelli, R. Siri, and V. Micenmacher. Phys. Rev. E, 
76:037202, 2007. 

8. A. Romanelli. Phys. Rev. A, 76:054306, 2007. 

9. Y. Yin, D. E. Katsanos, and S. N. Evangelou. Phys. Rev. A, 
77(2):022302, 2008. 

10. M. Stefandk, T. Kiss, I. Jex, and B. Mohring. /. Phys. A: Math. 
Gen., 39:14965, 2006. 

11. T. Kiss, I. Jex, and M. Stefanak. New Journal of Physics, 
11:043027, 2009. 

12. M. Stefaiiak, I. Jex, and T. Kiss. Phys. Rev. Lett., 100:020501, 
2008. 

13. T. Kiss, M. StefaMk, and I. Jex. Phys. Rev. A, 78:032306, 2008. 

14. N. Konno. Quan. Inf Proc, 9(3):405-418, 2010. 

15. N. Konno. Math. Struct. Camp. Sci, 20(6, Sp. Iss. SI):1091- 
1098, 2010. 

16. N. Konno. Quantum Information Processing, 8(5):387-399, 
2009. 

17. N. Konno. Quan. Inf & Comp., 2(Sp. Iss. SI):578-595, 2002. 

18. N. Konno. In H. Umeo, S. Morishita, K. Nishinari, T. Ko- 
matsuzaki, S. Bandini, editor. Cellular Automata, Proceedings, 
volume 5191 of Lectore Notes In Computer Science, pages 12- 
21, 2008. 

19. N. Konno. In M. Schurmann, U. Franz, editor. Quantum Po- 
tential Theory, volume 1954 of Lectore Notes In Mathematics, 
pages 309+, 2008. 

20. N. Konno. Stochastic Models, 25(l):28-49, 2009. 

21. N. Konno. Phys. Rev. E, 72(2, Part 2), 2005. 

22. N. Konno. Flue. And Noise Lett., 4:L529-L537, 2005. 

23. W. Diir, R. Raussendorf, V. M. Kendon, and H. J. Briegel. 
Phys. Rev. A, 66(5), 2002. 

24. B. C. Travaglione and G. J. Milburn. Phys. Rev. A, 65(3, Part 
A), 2002. 

25. 

26. M. Karski, L. Foerster, J.-M. Choi, A. Steffen, W. Alt, 
D. Meschede, and A. Widera. Science, 325(5937):174-177, 
2009. 

27. H. Schmitz, R. Matjeschk, Ch. Schneider, J. Glueckert, 
M. Enderlein, T. Huber, and T. Schaetz. Phys. Rev. Lett., 
103(9), 2009. 

28. F. Zaehringer, G. Kirchmair, R. Gerritsma, E. Solano, 
R. Blatt, and C. F. Roos. Phys. Rev. Lett., 104(10), 2010. 

29. A. Schreiber, K.N. Cassemiro, V. Potocek, A. Gabris, P. J. 
Mosley, E. Andersson, I. Jex, and Ch. Silberhorn. Phys. Rev. 
Lett., 104(5), 2010. 

30. M. El Ghafar, P. Torma, V. Savichev, E. Mayr, A. Zeiler, and 
W. R Schleich. Phys. Rev. Lett, 78(22):4181^184, 1997. 

31. S. W. Kim and H.-W. Lee. Phys. Rev. E, 61(5):5124r-5128, 
2000. 

32. R Ribeiro, R Milman, and R. Mosseri. Phys. Rev. Lett., 93(19), 
2004. 



33. G. Leung, P Knott, J. Bailey, and V. Kendon. New J. Phys., 
12(12):123018, 2010. 

34. H. Lavicka. Simulations of Agents on Social Networks. LAP 
Lambert Academic Publishing, 2010. 

35. J. Blank, P. Exner, and M. Havlicek. Linear Operators in 
Quantum Physics (in Czech). Carolinum Prague, 1993. 

36. M. Berry, I. Marzoli, and W. P. Schleich. Phys. World, 
14(6):39-44, 2001. 

37. R W. Anderson. Phys. Rev., 109(5):1492-1505, 1958. 

38. C. Tsallis. /. Stat. Phys., 52:479, 1988. 

39. C. Tsallis and E. Brigatti. Cont. Mech. and Ther., 16(3):223- 
235, 2004. 

40. C. Tsallis. Eur. Phys. f. A, 40(3):257-266, 2009. 

41. C. Tsallis. /. of Comp. and App. Math., 227(1, Sp. Iss. SI):51-58, 
2009. 

42. C. Tsallis. In W. Hillebrandt and F. Kupka, editor. Interdisci- 
plinary Aspects of Turbulence, volume 756, pages 21-48, 2009. 

43. D. Challet and Y.-C. Zhang. Physica A, 246(407), 1997. 



